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Abstract. We derive the equations of second order dissipative fluid dynamics from the 
relativistic Boltzmann equation following the method of W. Israel and J. M. Stewart [l]. 
We present a frame independent calculation of all first- and second-order terms and their 
coefficients using a linearised collision integral. Therefore, we restore all terms that were 
previously neglected in the original papers of W. Israel and J. M. Stewart. 

1 Introduction 

Relativistic fluid dynamics has been applied successfully to describe the dynamics of the hot and dense 
matter created in relativistic heavy-ion collisions at the Relativistic Heavy Ion Collider (RHIC) Q. It 
is also expected to play an important role in understanding of future experiments at the Large Hadron 
Collider (LHC) and at the Facility for Antiproton and Ion Research (FAIR). 

On the other hand, the theoretical foundations of relativistic dissipative fluid dynamics are not fully 
established yet. Relativistic kinetic theory of dilute gases provides a framework for relativistic fluid 
dynamics which can be derived systematically from the Boltzmann equation. However, this procedure 
is not unique and has been subject of many past and recent studies. 

In this paper we present a short but self-contained derivation of dissipative fluid dynamics from 
kinetic theory following the widely used approach by Israel and Stewart (IS) which is a generalization 
of Grad's method of moments |3| to relativistic systems. We show that the resulting equations contain 
several new second order terms and coefficients |4| which were absent in the treatment by IS and 
others in the past. 

Notation: We define the space-time coordinates choosing c = 1, and the co variant and contravari- 
ant four- vectors, - g^yx' = (t, -x, -y, -z) and x'' - g^^'x^ = (f, x, y, z), where g'^'' = - 
diag(l, -1, -I, -1) is the metric of flat space-time. The normalized hydrodynamic four-velocity of 
matter is u'^, such that m''m^ - 1. The transverse projection operator zl'^^ = g'"' - u^u^, is used to de- 
compose four- vectors or tensors into parts parallel and orthogonal to m'^. The transverse projection of 
four-vectors is denoted by - AyA^. The transverse and traceless projection of second-rank tensors 
is defined as A<'^''> = /1(^^"^ = [5 [a'^^A^^ + A^A''^) - ^zl^'M^^] A"^. The gradient, = d/dxf, of an 
arbitrary tensor can be decomposed as, cJ^A''' - ''" = UaDA^'^- ^^" + Vq, A'' ' where the comoving time- 
derivative, D = uf'di^,, is also denoted by an over-dot, = DA^, and V(j, = A^djj is the gradient operator. 
The symmetric and antisymmetric parts of a second rank tensors are denoted by A^'^^ = (A''^ + A'''^) /2 
and A^'^^ = (A''^ - A^'') /2, respectively. Using the above notations the relativistic generalization of the 
Cauchy-Stokes decomposition is dftUy - u^iiy + ^OAf^y + cr^v + co^y, where we defined the expansion 
scalar = V^m'^, the shear tensor cr''^ = V^u^\ and the vorticity tensor af^ = W^u^K 
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2 Fluid dynamics from l<inetic thieory 

In relativistic kinetic theory of single component gases, particles with mass m and four-momenta 
- (p^,P) where = -y/p^ + m?, at space-time coordinate - (t,x) are characterized by the 
invariant single particle distribution function / - fit, x, p). Assuming that there are no external forces, 
the space-time evolution of / in a dilute gas with binary collisions is given by the relativistic Boltzmann 
transport equation Q, 

7/5,/ = C[/], (1) 
where C [/] is the collision integral. The explicit form of the collision integral is 

cm = w,.,.„,/,,/,)|r/;//, -//,/■/;] , (2) 

where du = gcPpI [/7''(27rS)^j; g is the number of internal degrees of freedom and the transition rate 

is proportional to the cross section W{p,p\,p' ,p\) ~ cr. Here f = 1 - af denotes the correction due 
to boson and fermion statistics with a = -1 and a = +1, while a - Q for classical gases. 

For later convenience let us separate the four-momenta of particles into two parts using an arbitrary 
fluid dynamic flow velocity u''. Thus p^ - Eu^ + where E - p^u^ is the energy of the particle in 
the Local Rest Frame (LRF), where M^j^p = (1, 0, 0, 0), and = -^'^"pv is the LRF momentum. The 
macroscopic fields such as the particle four-flow and energy-momentum tensor are defined as the first 
and second moment of the single particle distribution function, 

N''(t, x) = Qy^) = {E)u'' + ip^i"^) , (3) 
T^'iU X) ^ ip^p') = {E^)u^^u' + jA^'i^^PaPp) + {Ep<''^X + {Ep^'^)u'' + {p^p'^) , (4) 

where the brackets, (...) = j doj . . . /, denote the integral of / in momentum space. Similarly, we can 
define higher moments or fluxes, these are F''^--''" = (p''^ . . . p^" ), and they obey the recurrence relation 
F^''■■■^^•'g^^^_^^^^ - ,7^2^/i,.../i„_2^ where p^p^ - was used. 

The conservation of particle number or charge and of energy and momentum in individual colli- 
sions lead to the equations of relativistic fluid dynamics ||5l . 

d^N'' = Jt/wC[/]=0, d^r'"" = ^ dcjp^C[f]^0, (5) 

while the divergence of higher moments leads to the balance of fluxes, 

g^pt>,...M„A = J j^pMi jjf„c [/] = pf^-f'.' , (6) 

where the pf denotes the production term for the n-th moment of the Boltzmann equation. Simi- 
larly to F'^i -'^", the production term also satisfies the recurrence relation, |^^^ = ,^j2p/ii.../i„-2 

The above equations are valid for an arbitrary solution / of the Boltzmann equation. However, in 
case the microscopic time and length scales are much shorter than the macroscopic ones, the system 
relaxes to the so-called local equilibrium state. In local equilibrium the single particle distribution 
function has a specific isotropic form resulting from a locally maximal entropy with vanishing entropy 
production, 

fo = [exp(-ao + /3oE) + a]-\ (7) 

where ySo - l/T is the inverse temperature and oq = p/T is the ratio of chemical potential /i and 
temperature, both in units of A;^ = 1. Using the local equiUbrium distribution function in Eqs. ( 1 3 141 ) we 
define, 

«o = <£>o, eo^{E\, po = -^<^"W/i>o, (8) 
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where the brackets, (. . .)() = j du . . . /o, and no - N^u^ is the particle number density, eo - T^^u^Uy is 
the energy density and po — -T'^A^^yl3 is the isotropic pressure, in equilibrium. Here we introduced 
the particle four-current N^ifo) = and energy-momentum tensor r'^''(/o) = ^y". Thus the local 
equilibrium approximation defines a perfect fluid with only five unknown fields m'', hq and eo, such 
that 

A'^ = nou" , (9) 
T^'' = eou^'u" - poA^'\ (10) 

where the pressure po is defined through an Equation of State (EOS). The temperature, chemical 
potential and entropy are also defined by the laws of equilibrium thermodynamics. 

If the system is out-of equilibrium the distribution function changes. Let us denote this diff'erence 
by, 5/ - f - fo, therefore the averages in Eqs. (1314b define 14 new fields on the l.h.s.. 



n = (E) = (£)o + {E}s ^no + 6n, (1 1) 

e = {E^}^{E\ + {E\,^eo + 6e, (12) 

p = = -li'^^^PaPdh - \{A''^PaPp)s ^PO+n, (13) 

y^^ = (14) 

W = = <£/?<^>)5 , (15) 

;i^'' = (pV>) = <P<V),, (16) 



where = / do) . . .6f, and n - N'^u^ - hq + 5n is the particle number density, e - T^"'u^Uy — 

eo + 6e is the energy density and p - po + TI - -T^^Af,y/3 is the isotropic pressure decomposed 
into equilibrium pressure and bulk viscous pressure parts. Here, V = A'^"Na is the particle diff'usion 
current with V^m^ - 0, W - A^^"Tapul^ is the energy-momentum diffusion current with Wm^ - 0, 
while the orthogonal and traceless part of T^^ defines the stress tensor, n^^' - t''^^'\ where tt'^'m^ - 
and n^'^g^jy = 0. Now, we can write the fundamental fluid dynamical quantities as 

s N^^ + SNf" = {no + Sn)u'' + V , (17) 
T"" = r^" + ST"" = (eo + 6e)ui'u' - (po + mA"" + IW^^u"^ +71^" . (18) 

At this point the flow of matter m'^ is usually specified using the definition of Eckart ^ or of Lan- 
dau and Lifshitz |7 |. Eckart chose the flow of conserved particles (if there are any) to define - 
/yjNf^Nfj, while Landau and Lifshitz chose the flow of energy-momentum = T^^Uyl^Ti'"Ua T^p . 
This implies that either the particle diffusion current or the energy-momentum diffusion current van- 
ishes, y = or W - 0. These physically different choices are related to each other and define the 
heat-flow, 

= ^y/^ _ v^'{e + p)/n - Vieo + po)/no + O2 • (19) 
In the following we will not restrict the calculations by fixing a frame of reference. 

3 The method of moments 

One of the methods successfully used to derive the equations of dissipative fluid dynamics was intro- 
duced by H. Grad B3I8I and later generalized to relativistic systems by Israel and Stewart 01I9I1OI1 II . 
This method uses the conservation laws together with the balance equation for the third moment and 
its production term to obtain the equations of motion for the 14-fields. 

Let us specify the distribution function as f(y) = (e"*' + a)"', where the equilibrium distribution 
function is recovered when y - ijQ = -/3oE. If the system is not too far away from local equilibrium 
the argument of the equilibrium distribution function changes hy Sy - y - yo 1 , therefore 



f(y) = fiyo) + f'(yo)Sy + f"(yo)6y^/2 + ... + 0„, 



(20) 
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where f'iyo) - /o/o- The method of IS specifies the deviation from equilibrium in terms of a polyno- 
mial expansion in the momentum four-vector p'^ up to quadratic order, 



6y = a-i//5,+ifp'' 
where /5 - P^Uf, and w 



uf^'/lfiv are scalars, if = jS^''^ and uf 



orthogonal to m'^, and u;^"^ is the traceless and orthogonal part of u/"'. 
Substituting the previous two equations into Eqs. ( I3I4I) we get 

N" = [ho + aJu) -/3J20 + w(Jio - J31)] -J2iif + IJiiuf , 



(21) 

^^l^up are four-vectors 



T"' = [/20 + a/zo -Phi) + w(J4Q - J41)] u^u' - 
- 2 (731 - 2741 w^^) m"^ + 2742if<''''^ . 



-h\ - aJii +pJn - If 1 741 - --/42 



(22) 



(23) 



Here we introduced the auxiliary thermodynamic integrals, I„q(ao,Po) and 7„^(q'o,j6o) 



where n and q < [n/2] are integers and {2q +1)11 denotes the double factorial. For a classical gas 
Jnqia = 0) = Ing. These two integrals are related by y6o7„, = -/«-i,^-i +(n- 2q')/„_i,^, dl„g = J„gdao - 
Jn+i,qdPu, andPodJnk - [-J„-i,q-i +(n- 2g)7„_i,5 f^ao+[-/«,?-i - (n - 2^ -t- 1)7„,] dPo where d stands 
for both the proper time and spatial derivatives IfTSI. 

Furthermore, comparing Eqs. ( fTTlfTSl l with Eqs. (l22ll23T l we can express the fluid dynamical fields 
in terms of the thermodynamic integrals and the parameters of the distribution function. 



«() = ^10 , eo - ho , Po - 
6n - aJio -pJio + wiho - hi) , 



-hi 



5e - aho -PJ30 + wUw - hi). 



U = -Q'721 -l-y673l - w(74l - -742 



y = -J2iif + 2hiuf .. 

71^" = 2742lf^''^ . 



W = -hiv^ + 2hiuf : 



(25) 
(26) 

(27) 

(28) 
(29) 



The above relations can be inverted to extract the parameters of the non-equilibrium distribution func- 
tion in terms of the 14 fluid dynamical fields, n - no + 5n, e = eo + Se, V , W, tt'"' and p = po + U. 
Without loss of generality the scalar quantities can be inferred from Eqs. ( I25I26I I. however the current 
method does not provide enough independent equations for closure. Therefore, one usually employs 
a so-called matching or fitting condition where the non-equilibrium particle and energy densities are 
assumed to be unchanged at least up to first order in deviations from equilibrium. Hence, following IS 
we fix 6n - 6e - 0. This also means that the EOS is given by, po - po{eo, no), and the temperature and 
chemical potential are fixed by the equilibrium state. Alternatively one could use different matching 
conditions |13| or even different thermodynamical theories, see for example Refs. lil4.15.16J . Thus 
from the above equations and constraints we get 



if 



w = -^nwhl , 



where we introduced 

^na - ^nm [m^ - 4 (73o73i - 72o74i) /D20] , 

-3D20/7 



^n/} - ^nm [4(7io74i - 72o73i) /D20] 



4 [372i(73o73i - 72o74i) -I- 373i(7io74i - 72o73i) - 5742D20] 
Sv„ = 741/D31 , Sw„ = 731/D31 , Svtt, = 73i/(2D3i) , Bww 

C„w - (2742) ' , 



/2l/(2£»3l), 



(30) 
(31) 
(32) 

(33) 

(34) 

(35) 
(36) 
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and D„q — Jn-\,qJn+\,q ~ Jnq- 

Here we note that using the definition of heat-flow together with the definition of enthalpy per 
particle /i = (eo+/'o)/no = /si /72i, one can equivalently choose (/' = -S'y^V^ +S'gi,q^ and uf -S'q^,q^, 
where S'y^ = {JnY^, = J31/D31 = Swd and = J21 /(2D31) = Sww, however we prefer to write 
the equations in terms of V and W. 

3.1 The balance equations from the third moment 

The equations of motion for the nine dissipative fields can be calculated from the third or higher 
moment of the Boltzmann equation. Therefore using Eq. (|6]l for n > 2, the equations for bulk viscous 
pressure, heat and diff'usion currents, and stress tensor can be formally written as 

M^, . . . u^,d,F'''-^"' = M^, . . . M^, P"' , (37) 
zj;, u,,_... u^,d,F'''-^'-' = zi;, u^,_... P'^-"" , (38) 
<f^,M^, . . . u,„d,F'"-^'"' = zf;f^,«^3 . . . P'^-"" . (39) 

Here we follow the method of IS and augment the equations of fluid dynamics Eqs. Q which are 
Eqs. (|6]l for n = and n - I with the next equation for n - 2 that is BaF^^"^ - P^^ . Using the 
14-moment approximation, F^^^'^ and its production term P''" become functions of the nine dissipative 
fields. This means that one discards irreducible tensors of rank higher than two, which usually do not 
appear in fluid dynamics. Of course one could obtain the equations of motion from higher moments of 
the Boltzmann equation, however, we note that the resulting equations of motion would be formally 
identical but yield different transport coefficients as shown in Ref. [ 17 |. 

In the following we explicitly calculate the equations of motion for the dissipative fields as done by 
Israel and Stewart, but here we restore all the terms that they neglected. Hence, using Eqs. ( 13013 1I32| | 
the third moment leads to 

F^'"^ = iifp^p^) = (/30 + iI^aH) uI'u'u^ + 3 (/31 - 1A4/7/3) m'^zI"'*' + 2>ilf'^ W^^ u'' u^^ + 2>^\v'^ u" u^'' 

+ 3iff'^W^''/!''^> + 3-^2 V'^zf'"'^ + SilfW^u^'^ , (40) 

where we introduced the following variables 

-/3{)(-/30-/3I - -/2O-/41) + -/40(-/l0-/41 - -/2O-/31) - -/sii'ao 

</'4 , (41) 

72l(-/30-/31 - J20J41) + J3l(JloJ4l - J20J31) - 5J42D20/3 

= (-^21/51 - /31/41) /D31 , = -</5 , (42) 

^J' = -D41 /D31 , (AI = mV5 - -Ar/5 , (43) 

iA3 = ^52/^42 . (44) 

The equations of motion follow from the different projections of the balance equation of the third 
moment. Thus Eq. (l37l i leads the equation for the bulk viscous pressure 77, 

u^UyP"" = 730 + </'477 + 1A477 + ^730 - 2731 + ^<A4^j 6 

+ ^^(lJ/Yw^') - 2iAf IV^M^ + df.iiJ/^V") - lif/^Vf^u" - lij/^n^'df^Uy . (45) 
The vector equation Eq. ( 1381 ) for W and V is 

Mv^r^p"' = (730 - 273, + |«A4^) + W [^r + ((Af - i/^Y)e] + [.Af + - 



A>1 



1 

(l/^r - 'A2"')^''5vM'' + - lljlWdyUl" + (dyllJ3 - l//3Uy)7T'"' + Ma^v^"" • (46) 
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The equation for the stress tensor ;r'"' can be calculated from Eq. ( l39l l, 

P'^"^ = 2 J/31 - ^ifr4nj d'^^'u"^ + 2(.Ar - </'r)M^''W"> + 2((/rJ' - 1^1)11^^^ + 2 (w<''5''Vf + lA^^^'W^^) 

+ 2 (y<''5''Vr + Vid^'V^''^) + hTT^-" + {tt^^"^ + TT^-'e) + li/.y^dAu"^ . (47) 

Next we evaluate the production term using the rh.s. of the Boltzmann equation. Using the distri- 
bution function from Eq. ( l20l l in Eq. (|2]l, the linearized production term is given by 

Pfv ^ 1 J dcodLJidco'dco\ W(/?,/?i,p',/?'i)/o7i;o/o/i,oe-(^'o+»;,o+!'o+</i,o)/2 

X [p"'//^ + p'I'pf - p"i/ - p"y^^ = C^^'"^w,p , (48) 

where the collision tensor, C^™fi - Qiiiv){ap)^ jg symmetric upon the interchange of two incoming or 
outgoing particles. Furthermore, the collision tensor is traceless C^^ga/i = C^^^l^g^tv - 0, and obeys 
time-reversal symmetry, C^™P - Using these properties we decompose the collision tensor as. 



+ ABQU^''A''^^''uf^ + ^Af'^^A'^^" , (49) 



where Aq = C^^^UfjUyUaUp, Bq = C^^"^U(^Ay)(aUi3)l2>, Co = C^^^Af^^a^^^y. Therefore using the above 
equations we easily get the l.h.s. of the balance equations. 



Uf.UyP^"' = 4Aqw/3 = CnU , 
UyA'^^P^" = IBauf = CvV + CwW" , 

where Cn = -(4AoJ?l;7^/3), Cv = -2BqSvw, Cw = 2BoSw., and C„ = CoC™/5. 



(50) 
(51) 
(52) 



3.2 The relaxation equations 

The relaxation equations follow from the balance equations ( l45l l46l l47T i and the linearized collision 
integral, Eqs. ( 150151152b . Here we write the equations in a frame independent form since it is easy to 
re-write them in the Eckart frame, where V = and q'' - W, or in the Landau and Lifshitz frame, 
where W = and then q^' = -hV. 

The relaxation equation for bulk viscosity is 



where we introduced G„q - JnoJqO - Jn-\,oJq+\,o and 



-1 



■ (no£>30 + (eo + Po)G23 - PqJai) , 



CnD2o 

Pn = = 'A4£>2()/ (no£>30 + (eo + ^0)^23 -PqJa\) , 



Tn ^ - 



-1 
C'n 
-1 



G23 



5A) £» 



20 



-1 



-1 

Dm 



2^1 +13 



(53) 

(54) 
(55) 

(56) 
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The relaxation equation for the flow of energy-momentum and conserved charge is 



W^' - hV^ = -K 



- ^SiwWO + h^divVe - AwwWyO^" + hAvvVyO^" , 

PO PO 



where we used that Cy = -hCw and defined, 



¥l D31 



Tw 



Cw J31 Cw 

/3w =fioTw/K = -UY + 



eo+Po 



Cv hCw 



-1/5, ySo#4 j6o^41 

^ \3 3 dfio eo + po; 



Aww 



Cw 

Cw 



- 1 \ 



-1 /1A4 _^ PqJai 
3 eo + Pa 



-1 



Iqn ^ —W3 + 



c 



1 

Cw 



eo + Po 



-1 ld>//4 _|_^-l#4l 

3Cw\dao d^o)' 



-1 

Cw 



eo+Po 



^qn - ^ 
_ -1 

Ayv - — 



1 



dao dPo I 
5 5 



5 

- + 



4A)^41 



6iv = Pv 



_i_ 5 _ 

3 "T 



3 3(eo + /?o)< 
The equation of motion for the shear stress tensor is 



(60) 

(61) 
(62) 
(63) 
(64) 
(65) 
(66) 
(67) 



where 



T] ■■ 
Ann ■ 



Cn' 

^4 

1 



Tn = - 



Cn' 



-1 (difr^ 



je. = Tj(277) = -^3/(2/31), 

w 



Cn I dao 



2 , f^-l9^2 



TnW - -pr 

^n 

W 



I 2 

' m 



-1 



dao dPo 



T^5^^-^°^l' 



, ^2 , 

^n 



(68) 

(69) 
(70) 
(71) 
(72) 



Note that we made extensive use of the following identities, V^jSq = h ^Wao - Pqu^ and V^if/ = 

(^.^h-^§^)'^^ao-Pogu^. 

The above relaxation equations are usually written and solved in a form which is given dividing 
the relaxation equations by their respective relaxation times. This leads to coefficients which do not 
depend on the cross section. Here we list them in the ultrarelativistic limit, where m/T 0, e = 3p 
and 77 = 0, in both the Eckart and Landau and Lifshitz frames. 



Pw -Pv ^ — — - , Pww ^ 7 ' Pvv - ~ * 7 ' 



Tw 



(73) 
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and 



^ V , 3»o , V 3«o 2 



y3;r ^ T— , y3;rW = — ' i^fi „ ' PnV^ ^ -T;;^— ' ^2 = ^l^Pn , (76) 

4po T„ llfiopo T„ 3y6g79o 3 



4 Conclusions 

In this work we derived the equations of second order dissipative fluid dynamics with all first- and 
second-order terms from the Boltzmann equation using the 14-moment method with a linearized col- 
lision integral. We also expressed all coefficients multiplying the second order terms independent of 
the choice of frame and showed that some of these coefficients are different in different frames. 
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